The period of a classical oscillator 
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We develop a simple method to obtain approximate analytical expressions for the period of a 
particle moving in a given potential. The method is inspired to the Linear Delta Expansion (LDE) 
and it is applied to a large class of potentials. Precise formulas for the period are obtained. 
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In this letter we consider the problem of calculating 
the period of a unit mass moving in a potential V(x). 
Although it is possible to solve this problem analytically 
only in a few cases, depending upon the form of the po- 
tential, several methods to find approximate results have 
been devised in the past. Many of the techniques that are 
used to solve this kind of problems are based on a per- 
turbative expansion in some small parameter that ap- 
pears in the equations of motion. This is the case of 
the Lindstedt-Poincare method and of the multiple-scale 
method. Unfortunately the validity of these approaches 
is restricted to the domain of weak couplings and the se- 
ries obtained rapidly diverge when larger couplings are 
considered. Recently, one of the authors and collabora- 
tors devised a non-perturbative version of the Lindstedt- 
Poincare method, based on the ideas of the Linear Delta 
Expansion (LDE)0|, which allows to obtain verv accu- 
rate results in a wide class of non linear problems 
In this letter we propose a different method, also inspired 
by the LDE, whose application is much simpler and for 
which convergence to the exact result can be proven. 

Let us describe the method in detail. We consider a 
unit mass moving in a potential V(x). The total energy 
E = s y+V(x) is conserved during the motion. The exact 
period of the oscillations will be given by: 
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where x± are the inversion points, obtained by solving 
the equation E = V(x±). 

Only in few cases it is possible to evaluate the integral 
JIJ analytically. In the spirit of the Linear Delta Expan- 
sion (LDE) we interpolate the full potential V(x) with a 
solvable one Vo(x) 1 defined as Vg{x) = Vq(x) + 6 (V(x) — 
V (x)). 

We want to perform this interpolation without moving 
the inversion points; for this reason we ask that x± be 
the inversion points also of the potential Vq(x). As a 
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V2 

V£o-v (zf 



result, the energy Eq that the particle would possess if it 
was moving only in the potential Vq(x) will be given by 
E Q = V (x±). 

We are now in position to rewrite eq. Q as 



x + 
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yjE - V (x) + S[E -Eq- V(x) + V (x)} 



dx. 
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We notice that Eq. @ reduces to Eq. QJ for 6 = 1; 
for S = this formula yields the period of oscillation 
between the points x± in the potential Vq(x). We will 
treat the term proportional to d as a perturbation and 
expand in powers of <5. Since Vq(x) depends upon one or 
more arbitrary parameters (which we will indicate with 
A) a residual dependence upon these parameters shows up 
in the period when the expansion is carried out to a finite 
order. In order to eliminate such unnatural dependence 
we impose the Principle of Minimal Sensitivity (PMS) 
by requiring that 
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Finally, we can write explicitly the period by perform- 
ing an expansion in 5 and obtain 



T S = 



V2 



dx 



x _ v/£ - V Q (x) ^1 + 6 A(x) 
x _ y/Eo - V (x) 



where 
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E-E Q -V(x) + V {x) 



E -V (x) 



Tx can be written as 
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provided that the series in Eq. J3J converges uniformly, 
which is the case if |A(x)| < 1 for every x, X- < x < X+. 

We now apply Eq. JSJ to a few cases. We start with 
the Duffing oscillator, which corresponds to the potential 
V(x) = f) x 2 + j x 4 . Although the period of the Duffing 
oscillator can be calculated explicitly in terms of elliptic 
functions, we use this example to illustrate our method 
and to prove its efficiency. 

We choose the interpolating potential to be Vq(x) = 
1+ „ x x 2 and obtain 



A(x) 



1 + A 2 
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Hence the series in Eq. converges to the exact pe- 

since |A(x) 
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< 1 



riod for A > Ao = y ^j- y 1 
uniformly for such values of A and |a;| < A. 

The period of the Duffing oscillator calculated to first 
order using © is then 
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By setting 8=1 and applying the PMS we obtain the 
optimal value of A, Xpms = "*2 4i wmcn remarkably co- 
incides with the one obtained in |2j by using the LPLDE 
method to third order. The period corresponding to the 
optimal A is 
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which provides an error less than 2.2% to the exact period 
for any value of p, and A. 

It is useful to compare this result with the one in 0, 
which differs from our result only for a numerical factor 
under the square root. The result of Q yields a much 
larger error (16%) over the period. This is expected, 
since our first order formula complies with the PMS. 

Let us now come to the issue of convergence: since 
Xpms > Ao we can write 



Tt = 



(G) 



where 
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FIG. 1: Error over the period (in absolute value), defined 
x 100, for A = 10 and /j, = 1 as a 



function of the order. The th ree sets are obtained by using 
the optimal value Xpms = \/3 fj>A 2 /2 (plus), a value A = 
0.9 Xpms (triangle) and A = 1.1 Xpms (square). 



exact result, which is precisely what we observe in Fig- 



x 100 



ure^ where we plot the error S 
for three different values of the parameter A as a function 
of the order in the expansion. T exact is the exact period 
of the Duffing oscillator. Corresponding to the optimal 
value of the parameter, Xpms = y/3 yuA 2 /2, the rate of 
convergence is maximal. 

We now consider the general anharmonic potential 2 
V(x) = h x 2 -f- p Tny- and obtain 



A(x) 



1 + A 2 
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To first order, Eq. (0 gives 
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in which case the optimal value of A is given by 
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This time, these values of A do not coincide with the 
ones obtained with the LPLDE method to third order 
for N — 3 and TV = 4 Q. With such A one obtains the 
expression 

2- 
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2pA^ N ^T{N + 1/2) 
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2-Fi is the hypergeometric function. Since Eq. Q is es- 



2 The Duffing oscillator considered in the previous example corre- 
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We have tested Eq. JHJl for moderate values of p and A 
and seen that it provides a very good approximation to 
the exact period, even when the anharmonicity exponent 
N gets very large. 

We consider now the nonlinear pendulum, whose po- 
tential is given by V(9) = 1 — cos 6*. By choosing the 



interpolating potential to be Vq(9) 



1 + X±n2 



A(9) 



2 cos — cos ( 
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9 we obtain 
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where O is the amplitude of the oscillations. To first 
order our formula yields 
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where J\ is the Bessel function of the first kind of order 
1. The optimal value of A in this case is given by 



2Ji(6) 
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and the period to first order is then 
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Eq. (0 provides an excellent approximation to the ex- 
act period over a wide range of amplitudes. 

We now apply our expansion to two problems in Gen- 
eral Relativity: the calculation of the deflection of the 
light by the Sun and the calculation of the precession of 
a planet orbiting around the Sun. We use the notation 
of Weinberg Q : 



B(r) = A -1 (r) = 1 - 



2GM 



The angle of deflection of the light by the Sun is given 
by the expression 



A0 



B(r) 



1 -1/2 



dr 
r 



where r$ is the closest approach. 

With the change of variable z = 1/r we obtain 



Acj) = 2 r\ 
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l/r 



dz 
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2GM{\ -r^z 2 ) 



which is exactly in the form required by our method. We 
introduce the potential Vq(z) = (tq — A 2 ) z 2 to obtain 



A(z) = 



r (l - r 2 z 2 ) - 2GM (1 - r$z A ) 



( J. _ .2^ 



By applying our method to first order we obtain 
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The optimal value of the parameter A obtained by the 
PMS is given by 



2GM 



The deflection angle corresponding to this value of A 
reads 



&<t>PMS = -7r 
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1 - 



8GM 
r 7T 



(11) 



The surface corresponding to the closest approach for 
which Acf> diverges is known as photon sphere and for the 
Schwartzchild metric takes the value ro = 3GM. It is 
remarkable that Eq. , despite its simplicity, is able to 
predict a slightly smaller photon sphere, corresponding to 
ro = 8GM/w. Clearly, this feature is completely missed 
in a pcrturbativc approach. 

In FigureElwe compare Eq. 1)11(1 with the exact numer- 
ical result, the post-post-Newtonian (PPN) result of || 
and with the asymptotic result for very small values of ro 
(close to the photon sphere). In Figurc|31we plot the error 
over the deflection angle obtained by using Eq. (|11|) or 
the PPN result in a range of r of practical interest. Our 
formula gives an error which is two orders of magnitude 
smaller than the one obtained with the PPN. 

Finally we consider the problem of calculating the pre- 
cession of the perihelion of a planet orbiting around the 
Sun. The angular precession is given by 



A9 



where E 
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yjA(r) dr 



E 
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J 2 B(r) J 2 r 2 
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B(r+) S(r_) / 

(flfej- s^)/(^-pr)- r ± are the shortest 
(perielia) and largest (afelia) distances from the sun. By 
the change of variable z = 1/r we can write Eq. 1)12(1 as 



A9 = -2 



1 



V( z + - z)(z- Z-) 
dz 



y/(l - 2GM(z + Z_ + Z+)) 
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FIG. 2: Deflection angle of light obtained assuming G/c 2 



FIG. 4: Precession of the orbit of a planet assuming the values 



7.425 x 10 m/kg and M = 1.97 x 10 30 kg as function of the M = 1.97 x 10 30 kg, G/c 2 = 7.425 x 1CT 30 m/kg and e 
closest approach tq. The solid line is the exact (numerical) 
result, the dashed line is obtained with Eg. <f 1 1 p . the dotted 
line is the post-post-Newtonian result of the dot-dashed 
line is the asymptotic result (ro — » oo). The vertical line 
marks the location of the photon sphere, where the deflection 
angle diverges. 



0.2506 (eccentricity). The scale of reference oo is taken to be 
the semimajor axis of Mercury's orbit (ao = 5.971 x 10 10 m). 
The solid line is the exact result, the dashed line is the result 
of Eq. 11131 and the dotted line is the leading term in the 
perturbative expansion. 




FIG. 3: Absolute value of the error over the deflection angle 
as a function of the closest approach in units of the sun radius 
(r S un = 6.95 x 10 8 m), assuming G/c 2 = 7.425 x 10~ 30 m/kg 
and M = 1.97 x 10 30 kg. The error is defined as A = 
A ^" pp r° x — A * . The solid line is the error obtained with 

&f exact 

eq. 1111 , while the dashed line is the error obtained using the 
PPN approximation of 0. The vertical line marks the solar 
radius. 



We write 



A6 S = -2 



1 



- V( z + - z )( z - z -) 

dz 



y/(l - A 2 + 5(-2GM(z + z_ + z+) + A 2 )) 



2tt. 



As usual we treat the term proportional to 5 as a per- 
turbation and expand to third order. The optimal value 
of A is obtained by using the PMS 



6GM 



and yields the precession 
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tt(3G 2 LM 2 



-4L 2 + 48GLM ~ U7G 2 M 2 )) 



4a(L-6GM) 2 dl- 



6 G M 
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where a is the semimajor axis of the ellipse, given by 
a = (r_ + r_|_)/2, and L is the semilatus rectum of the 
ellipse, given by 1/L = (l/r+ + l/r_)/2. 

In Figure 01 we plot the precession of the orbit cal- 
culated through the exact formula (solid line), through 
Eq. 1)1 3|) (dashed line) and through the leading order re- 
sult A6> = §1 ^ M - (dotted line) . Once again we find 
excellent agreement with the exact result. 

In conclusion, we have devised a method to calculate 
with high accuracy a certain class of integrals, which are 
very common in many physical problems. The conver- 
gence of our expansion to the exact result is easy to verify, 
as we have explicitly shown in one special case. More- 
over the lowest order results obtained by applying our 
method already provide an excellent agreement with the 
full exact results. Work is currently in progress to apply 
this technique to a wider class of problems. 
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